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^ Abstract 

Wc study the problem of computing geometric spanners for (addi- 
Q tively) weighted point sets. A weighted point set is a set of pairs {p, r) 

'— ' where p is a point in the plane and r is a real number. The distance 

between two points {pi,ri) and (pj,rj) is defined as \piPj\ — ri — rj. 
^ We show that in the case where all rj are positive numbers and 

^ \piPj\ > Vi + rj for all i,j (in which case the points can be seen as 

Q non-intersecting disks in the plane), a variant of the Yao graph is a 

(1 + e)-spanner that has a linear number of edges. We also show that 
T-H the Additively Weighted Delaunay graph (the face-dual of the Addi- 

^ tively Weighted Voronoi diagram) has constant spanning ratio. The 

^ straight line embedding of the Additively Weighted Delaunay graph 

^ may not be a plane graph. We show how to compute a plane embed- 

• rH ding that also has a constant spanning ratio. 

X 



*Research partially supported by NSERC, MRI, CFI, and MITACS. 



1 



1 Introduction 



Let G be a complete weighted graph where edges have positive weight. Given 
two vertices u,v of G, we denote by Sg{u,v) the length of a shortest path 
in G between u and v. A spanning subgraph if of G is a t-spanner of G if 
Sh{u,v) < tdciu.v) for all pair of vertices u and v. The smallest t having 
this property is called the spanning ratio of the graph H with respect to G. 
Thus, a graph with spanning ratio t approximates the (2) distances between 
the vertices of G within a factor of t. Let P be a set of n points in the plane. 
A geometric graph with vertex set P is an undirected graph whose edges are 
line segments that are weighted by their length. The problem of constructing 
t-spanners of geometric graphs with 0{n) edges for any given point set has 
been studied extensively; see the book by Narasimhan and Smid p5j for an 
overview. 

In this paper, we address the problem of computing geometric spanners 
with additive constraints on the points. More precisely, we define a weighted 
point set as a set of pairs (p, r) where p is a point in the plane and r is a real 
number. The distance between two points {pi,ri) and {pj,rj) is defined as 
\piPj\—ri — rj. The problem we address is to compute a spanner of a complete 
graph on a weighted point set. To the best of our knowledge, the problem 
of constructing a geometric spanner in this context has not been previously 
addressed. We show how the Yao graph can be adapted to compute a (l + e)- 
spanner in the case where all are positive real numbers and \piPj | > Tj + rj 
for all i,j (in which case the points can be seen as non-intersecting disks in 
the plane). In the same case, we also how the Additively Weighted Delaunay 
graph (the face-dual of the Additively Weighted Voronoi diagram) provides 
a plane spanner that has the same spanning ratio as the Delaunay graph of 
a set of points. 

1.1 Motivations 

It has been claimed (see |2l [30l EI]) that geometric spanners can be used 
to address the link selection problem in wireless networks. In most cases, 
however, two assumptions are made: 

1. nodes can be represented as points in the plane and 

2. the cost of routing a message is a function of the length of the links 
that are successively used. 
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However, these assumptions do not always hold. For example, the first as- 
sumption does not hold in the case of wide area mesh networks, where nodes 
are vast areas such as villages [28]. The second assumption does not take 
into account the fact that some nodes may have higher energy resources or 
introduce more delay than others. In such cases, an additional cost must 
be taken into account for each node. The study of spanners of additively 
weighted point sets is a first step in addressing some of these issues. 

1.2 Paper Organization 

The rest of the paper is divided as follows: In Section [2| we review related 
work. In Section [3} we give a formal definition of our problem and show that 
it is not solved by a straightforward extension of the Yao graph. However, 
in Section |4], we show that a minor adjustment to the Yao graph allows 
to compute a (1 + e)-spanner. In Section [sj we develop some tools used 
in Section |6] to show that the Additively Weighted Delaunay graph has a 
constant spanning ratio. We conclude in Section [7) 

2 Related Work 

Well known examples of geometric t-spanners include the Yao graph [53] . 
the 6'-graph p9], the Delaunay graph [I9], and the Well-Separated Pair De- 
composition (WSPD) [8j. Let 6 < 7r/4 be an angle such that 2-n/6 = k, 
where k is an integer. The Yao graph with angle 6 is defined as follows. For 
every point p, partition the plane into k cones Cp^i, . . . , Cp^k of angle 9 and 
apex p. Then, there is an oriented edge from p to g if and only if q is the 
closest point to p in some cone Cp^i. The Yao graph is sometimes confused 
with the 6'-graph, although they are different graphs. The first phase of the 
construction of the 6'-graph using k cones with angle 9 and apex p is identical 
to the construction of the Yao graph. This may be the root of the confu- 
sion. However, there is an oriented edge from p to g in the ^-graph if and 
only if q has the shortest projection on the bisector of the cone containing q. 
For Yao graphs [33j, the spanning ratio is at most l/(cos6' — sin 6*) provided 
that 9 < Tc/A, and for ^-graphs, the spanning ratio is at most 1/(1 — 2 sin |) 
provided that ^ < 7r/3 [29]. 

Given a set of points in the plane, there is an edge between p and q 
in the Delaunay graph if and only if there is an empty circle with p and 
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Figure 1: The Delaunay graph and its dual: the Voronoi diagram. 



q on its boundary [19j. If no four points are cocircular, then the Delau- 
nay graph is a uniquely defined near-triangulation. Otherwise, four or more 
co-circular points may create crossings. In that case, removing edges that 
cause crossings leads to a Delaunay triangulation. Since our results hold for 
any Delaunay triangulation, when we refer to the Delaunay triangulation in 
the case of co-circular points, we mean any Delaunay triangulation. Dobkin 
et al. pjj showed that the Delaunay triangulation has a spanning ratio of 
at most i^TT ^ 5.08. This result was improved by Keil and Gutwin [T9] . 
who showed that the spanning ratio of the Delaunay triangulation is at most 
27r/ (3 cos(7r/6)) ~ 2.42. Later, Bose et al. [7] showed that the Delaunay trian- 
gulation is also a strong t-spanner for the same constant t = 27i/ (3 cos(7r/6)). 
Although the exact spanning ratio of the Delaunay triangulation is unknown, 
it is conjectured that the spanning ratio is tt/2. For the remainder of this 
paper, we will refer to the spanning ratio of the Delaunay triangulation as 
the spanning ratio of the standard Delaunay triangulation and denote it as 
SP-DT. 

The Voronoi diagram of a finite set of points P is a partition of the 
plane into \P\ regions such that each region contains exactly those points 
having the same nearest neighbor in P. The points in P are also called sites. 
It is well known that the Voronoi diagram of a set of points is the face dual 
of the Delaunay graph of that set of points [10], i.e. two points have adjacent 
Voronoi regions if and only if they share an edge in the Delaunay graph (see 
Figure [T]) . 

Let s > be a real number. Two set of points A and B in M"' are well- 
separated with respect to s if there exists two dimensional balls Ca and Cb 
of same radius r respectively containing the bounding boxes of A and B such 
that the distance between Ca and is greater than or equal to s x r. The 
distance between Ca and is defined as the distance between their centers 
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minus 2 x r. A Well- Separated Pair Decomposition with separation ratio s of 
a set of points P [HI ES] is a set of unordered pairs {{^i, -Bi}, . . . , {Am, Bm}} 
of subsets of P that are well-separated with respect to s with the additional 
property that for every two points p,q E P there is exactly one pair {Ai, Bi} 
such that p E Ai and q G B^. Callahan and Kosaraju [H] showed that for 
s > 4, every point set admits a WSPD with separation ratio s of 0{n) size 
that can be computed in 0(n\ogn) time. Choosing one edge per pair allows 
to construct a t-spanner that has 0{n) size with t = (s + 4)/(s — 4). 

Our work falls in the context of computing spanners for geometric graphs 
other than the complete Euclidean graph. Typically, variations of the spanner 
problem arise by either changing the distance function or removing edges 
from the complete graph. For example, for a set P of points in the plane 
and a set C of non-intersecting line segments whose endpoints are in P, the 
visibility graph of P with respect to C is the geometric graph with vertex set 
P and there is an edge [pq) if and only if the segment pg is in C or it does 
not cross any segment in C (in that case, p and q are said to be visible). A 
spanner of the visibility graph should then approximate Euclidean distances 
for every pair of points that are visible from each other. The constrained 
Delaunay triangulation (a variation of the Delaunay triangulation) is a 2.42- 
spanner of the visibility graph [6l [T71 [22] . 

Unit disk graphs [HI HU] received a lot of attention from the wireless 
community. A unit disk graph is a graph whose nodes are points in the 
plane and edges join two points whose distance is at most one unit. It is 
well-known that intersecting a unit disk graph with the Delaunay or the Yao 
graph of the points provides a t-spanner of the unit disk graph [7j, where the 
constant t is the same as the one of the original graph. However, this simple 
strategy does not work with all spanners. In particular, it does not work 
with the ^-graph [5]. Unit disk graphs can be seen as intersection graphs 
of disks of same radius in the plane. The general problem of computing 
spanners for geometric intersection graphs has been studied by Purer and 
Kasiviswanathan [THj . 

Another graph that has been looked at is the complete k-partite Euclidean 
graph. In that case, points are assigned a unique color (which may be thought 
of as a positive integer) between 1 and k, and there is an edge between two 
points if and only if they are assigned different colors. Bose et al. |1] showed 
that the WSPD can be adapted to compute a t-spanner of that graph that 
has 0{n) edges for arbitrary values of t strictly greater than 5. 

Por spanners of arbitrary geometric graphs, much less is known. Althofer 
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et al. have shown that for any t > 1, every weighted graph G with n 
vertices contains a subgraph with 0(n-'^+^/*^*~^)) edges, which is a t-spanner 
of G. Observe that this result holds for any weighted graph; in particular, 
it is valid for any geometric graph. For geometric graphs, a lower bound 
was given by Gudmundsson and Smid [T3]: They proved that for every real 
number t with 1 < t < ^logra, there exists a geometric graph H with n 
vertices, such that every t-spanner of H contains ^2(72^+^/*) edges. Thus, if 
we are looking for spanners with 0{n) edges of arbitrary geometric graphs, 
then the best spanning ratio we can obtain is G(logn). 

In the literature, spanners that use a distance other than the Euclidean 
distance have also been proposed. For example, in a power spanner [21 dH 
EH inn], the distance used to measure the length of an edge is the square 
of the Euclidean distance between its two end points. This models the fact 
that in wireless networks, the amount of energy needed to send a packet is 
proportional to a power (not necessarily the square, however) of the Euclidean 
distance between the sender and receiver When reducing the latency 
is more important than reducing the amount of energy being used, a hop 
spanner p], which gives an equal weight to every edge, can be used. 

In this paper, the Additively Weighted Voronoi diagram (AW-Voronoi di- 
agram) is of particular interest. Its definition is similar to that of the (stan- 
dard) Voronoi diagram, except that each site Pi is assigned a weight which 
is a real number r^. Weights are used to define a weighted distance. More 
detail about how the weighted distance is used to define the AW-Voronoi 
diagram is given in Section [6j The Additively Weighted Delaunay graph 
(AW-Delaunay graph) is defined as the face-dual of the AW-Voronoi dia- 
gram. Properties of the AW-Voronoi diagram and its dual have been studied 
by Lee and Drysdale [22], who showed how to compute it in 0{n log^ n) time. 
Later on. Fortune [12] showed how to compute it in 0(n log n) time. The 
AW-Voronoi diagram may have empty cells. For this reason, one would hope 
that it is possible to design an algorithm whose running time gets better as 
the number of empty cells increases. Karavelas and Yvinec [TH] provided 
an 0{nT{h) + hlogh) time algorithm to compute the AW-Voronoi diagram 
where h is the number of non-empty cells and T{h) is the time to locate the 
nearest neighbor of a query point within a set of h points. Experimental 
results suggested an 0{n\ogh) behavior. In 3D, the complexity of the (Ad- 
ditively Weighted) Voronoi diagram is 0(?7.^) [21]. Aurenhammer [3] showed 
how to compute it in time 0(r;,^) using Power Voronoi diagrams. Will [32] 
gave an 0(n^ log n) time algorithm with experimental results suggesting an 



6 



0{n\og'^ n) time behavior in the expected case. Kim et al. [20] showed how 
to obtain a running time of 0{nm), where m is the number of edges. 

3 Definitions and Notation 

Definition 3.1 A set P = {(pi, ri), . . . , r„)} of ordered pairs, where 
each Pi is a point in the plane and each ri is a real number, is called a 
weighted point set. The notation pi E P means that there exists an ordered 
pair {pi,ri) such that {pi,ri) G P. The additive distance from a point p ^ P 
in the plane to a point Pi G P , noted d{p,pi), is defined as \ppi\ — ri, where 
\ppi\ is the Euclidean distance from p to pi. The additive distance between 
two points Pi,Pj G P, noted d{pi,pj), is defined as \piPj \ —ri — rj, where \piPj\ 
is the Euclidean distance from pi to pj . 

The problem we address in this paper is the following: 

Problem 3.2 Let P be a weighted point set and let K{P) be the complete 
weighted graph with vertex set P and edges weighted by the additive distance 
between their endpoints. Compute a t-spanner with 0{n) edges of K{P) for 
a fixed constant t > 1. 

Notice that in the case where all rj are positive numbers, the pairs {pi, r^) 
can be viewed as disks Di in the plane. If, for all i, j we also have d{pi, pj) > 0, 
then the disks are disjoint. In that case, the distance d{Di, Dj) = d{pi,pj) = 
\piPj \ — Tj — Tj is also equal to min{|gjgj| : qi E Di and qj G Dj}, where the 
notation G Di means \piqi\ < r^. To compute a spanner of an additively 
weighted point set is then equivalent to computing a spanner of a set of disks 
in the plane. From now to the end of this paper, it is assumed that 
all Tj are positive numbers and d{pi,pj) > for all i,j. If P is a set of 
disks in the plane, then a spanner of P is a spanner of the complete graph 
whose vertex set is V and whose edges {Di, Dj) are given weights equal to 
diDi,Dj). 

Notice also that the additive distance may not be a metric since the 
triangle inequality does not necessarily hold (see Figure [2]). Although this 
may seem counter-intuitive, this makes sense in some networks, since a direct 
communication is not always easier than routing through a common neighbor. 
For example, in wireless networks, the amount of energy that is needed to 
transmit a message is a power of the Euclidean distance between the sender 
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Figure 2: The additively weighted distance is not a metric. 

and the receiver. Therefore, using several small hops can be more energy 
efficient that a direct communication over one long-distance link. 




Figure 3: A straightforward generalization of the Yao graph. 

Figure |3] shows how the Yao graph can be generalized using the additive 
distance: every node keeps an outgoing edge with the closest disk that inter- 
sects each cone. However, this graph is not a spanner. Figure |4] shows how to 
construct an example with four disks that has an arbitrarily large spanning 
ratio. Nonetheless, in Section |4| we see that a minor adjustment to the Yao 
graph can be made in order to compute a (1 + e)-spanner of a set of disjoint 
disks that has 0(n) edges. 
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Figure 4: The straightforward generahzation of the Yao graph does not have 
constant spanning ratio. 




The Delaunay graph in the additively weighted setting is computable in 
time O(nlogn) [12J. To the best of our knowledge, its spanning properties 
have not been previously studied. In the two next sections, we show that 
it is a spanner and that its spanning ratio is SP-DT (i.e the same as that 
of the standard Delaunay graph). Finally, we show that although the addi- 
tively weighted Delaunay graph is not necessarily plane, it contains a plane 
subgraph that is a spanner with the same spanning ratio. 



4 The Additively Weighted Yao Graph 

As we saw in the previous section, a straightforward generalization of the Yao 
graph fails to provide a graph with bounded spanning ratio. In this section. 
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we show how a few subtle modifications to the construction, provide an 
approach to build a (l + e)-spanner. We define the modified Yao construction 
below. 

Definition 4.1 Let V be a finite set of disjoint disks and 9 < 0.228 he an 

angle such that 2ti/9 = k, where k is an integer. The Yao(6',P) graph is 
defined as follows. For every disk D = {p, r), partition the plane into k cones 
Cp^i, . . . , Cp^k of angle 9 and apex p. A disk blocks a cone Cp^i provided that 
the disk intersects both rays of Cp^i. Let F E T> he a disk different from D 
with center in Cpj. Add an edge from D to F in Yao{9, V) if and only if one 
of the two following conditions is met: 

1. among all blocking disks that have their center in Cpj, F is the one 
that is the closest to D; 

2. among all disks that have their center in Cpj and are at a distance of 
at least r to D, F is the one that is the closest to D. 

Notice that there are two main changes. Within each cone, we now add 
potentially two edges as opposed to only one edge in the case of unweighted 
points. Next, in the second condition to add an edge, we do not add an edge 
to the closest disk within a cone but to the closest disk whose distance is at 
least r from the disk centered at the apex with radius r. We now prove that 
these two modifications imply that the resulting graph is a (1 + e)-spanner. 

Lemma 4.2 Let pi,p2,P3 such that the angle Zp3pip2 = a < 9 < it /A and 
biPsI < Ipml- Then jpaPsI < Ipml - (cos6' - sin6')|piP3|. 

Proof: Let p'^ be the projection of ps on the line through pi and p2 (see 
Figure |5]). Then 

\P2P3\ < \P2P3\ + Ip'sPsI 

= \piP2\ - Ipip'sI + bsPsI 

= \PiP2 1 - \piP3 1 (cos a -sin a) 
< \PiP2 1 - \PiP3 1 (cos 9 - sin 9) 

□ 



Theorem 4.3 Let V be a finite set of disjoint disks and 9 < 0.228. Then 
Y{9,V) is a t-spanner of T> , where t = 1/ {cos29 — sm.29 — 2sin{9 /2)). 
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Figure 6: If D2 blocks the cone but the edge {Di,D2) is not in Yao{9,V), 
then there exists such that d{Di,D^) + d{D^, D2) < d{Di, 02). 

Proof: We proceed by induction on the rank of the weighted distances 
between the pairs of disks Di and D2. 

Base case: The disks Di and D2 form a closest pair. In that case, the 
edge iDi,D2) is in Yao{e,V). 

Induction case: Let Di = (pi,ri) and D2 = {p2,f2)- Without loss of 
generality, ri < r2. If the edge {Di, D2) is in Yao(6', V), then there is nothing 
to prove. Otherwise, there are two cases to consider depending on whether 
or not the shortest path from Di to D2 in the complete graph on V is the 
edge {Di,D2). If the shortest path is not the edge (1)1,1)2), then all edges 
on the shortest path must have length less than d{Di, D2). By applying the 
induction hypothesis on each of those edges, we conclude that the distance 
from Di to D2 in Yao(6', V) is at most t times the length of the shortest path 
Di to D2 in the complete graph on V, as required. 

We now consider the case when the edge {Di, D2): 

1. is not in Yao(6',P) and 

2. is the shortest path from Di to D2 in the complete graph. 

Observe that the conjunction of those two facts imply that the disk D2 
does not block the cone whose apex is pi and contains p2'- If -D2 was blocking 
the cone, then since (Z)i,Z)2) is not an edge in Ya.o{9,V), there must be a 
disk that is also blocking the cone and is closer to Di than D2. However, 
this implies that the shortest path from Di to D2 in the complete graph is 
not the edge (Di, D2) (see Figure |6]). 

The conjunction of the three following facts: 

1. ri < r2; 
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Figure 7: Illustration of the proof of Theorem 



4.3 



2. e < 0.228 < sin-^(l/3) and 

3. D2 does not block the cone, 

imply that d{Di, D2) > ri. Since {Di, D2) is not an edge, there another disk 
whose distance is at least r that is closer to Di. Let = (p3,r3) be the 
closest disk to Di such that ps is in the same 6'-cone with apex at pi as p2 and 
d{Di,D^) > ri. By definition, the edge {Di,D3) is in Yao(6',r'). Observe 
that d{D2,D^) < d{Di,D2). To see this, let a := d{Di,D2) - n. We have 
that 

d{D2, D3) <a + 4ri sin(^/2) < a + 4ri sin(0.114) < a + ri = d{Di, D2). 

Let p[ be the point of Di that is the closest to D^, p'[ be the point of Di 
that is the closest to D21 p'2 be the point of D2 that is the closest to Di, 
and P3 be the point of D^, that is the closest to Di (see Figure [?]). Notice 
that IP1P3I < IP1P2I that since d{Di,D2) > d{Di,Ds) > ri, then the 
angle /-p'^p'xP^j, is at most 1Q < tt/4. Therefore, we can apply Lemma 4.2 to 
conclude that 

IP2P3I < IP1P2I - (cos 2^ - sin2^)|p;p3|, 

which implies that 

d{D2,D3) < d{Di,D2) + \p[p"\ - {cos2e ~sm2e)d{Di,D3). 
Also, since \p\p'l\ < 2sin(0/2)ri < 2sm{e /2)d{Di, D3), we have 

d{D2,D3) < d{Di,D2) - (cos 2^ - sin 2^ - 2 sin(^/2))c/(L)i, D3). 
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Finally, since d{D2,D3) < d{Di,D2), the induction hypothesis tells us 

that Yao(^^, V) contains a path from D2 to D3 whose length is at most 
td{D2, D^). This means that the distance from Di to D2 in Yao(6',P) is 
at most 

d{Di,D3)+td{D2,Ds) < d{D,,Ds)+t{d{Di,D2)-^d{D^,D3))^td{Di,D2). 
The value 0.228 is an upper bound on the values of 9 such that t > 0. □ 

Corollary 4.4 For any e > and any set Vofn disjoint disks, it is possible 
to compute a (1 + e)- spanner ofV that has 0{n) edges. 

Proof: The bound on the number of edges comes from the fact that each 

cone contains at most two edges, and the stretch factor of 1 + e comes from 
the fact that lim l/(cos2^ - sin 2^ - 2sin(^/2)) = 1. □ 



5 Quotient Graphs and Quotient Spanners 

The main idea in the remainder of this paper is the following: we show how 
to compute a set of points from each such that the (standard) Delau- 
nay graph of those points is equivalent to the Additively Weighted Delaunay 
graph. By choosing the appropriate equivalence relation as well as the appro- 
priate point set, we can then show that the spanning ratio of the Additively 
Weighted Delaunay graph is bounded by the spanning ratio of the standard 
Delaunay graph. The reduction of one graph to another is done by means of 
a quotient: 

Definition 5.1 Let Pi and P2 be non-empty sets of points in the plane. The 
distance between Pi and P2, denoted by IP1P2I, is defined as the minimum 
\p\P2\ over all pairs of points such that pi e Pi and p2 G P2. 

Definition 5.2 Let G = {V,E) be a geometric graph and V be a partition of 
V. The quotient graph of G by V, denoted G/V, is the graph having V as 
vertices and there is an edge {U, W) (where U and W are in V) if and only 
if there exists an edge {u, w) G E with u & U and w e W. The weight of the 
edge {U,W) is equal to \UW\. 
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Del(P)/7' 



Figure 8: Illustration of Lemma 5.3 



If P is a (non-weighted) point set and P is a partition of P, then the 
notation P/V designates the quotient of the complete Euclidean graph on P 
by V. If 5 is a set of pairwise disjoint sets of points in the plane such that 
P C IJ5, then the notation P/S designates the quotient of the complete 
Euclidean graph on P by the partition of P induced by S. 

Lemma 5.3 Let G = {V, E) be a complete geometric graph, V be a partition 
of V and S be a t-spanner of G. Then S/V is a t-spanner of G/V. 

Proof: Let {U, W) be an edge of G/V and {u, w) be an edge of G such that 
\uw\ = \UW\. Since G is complete, the edge {u,w) is in G, and since S' is a 
t-spanner of G, there is a path = mi, . . . , in such that ui = u,Uk = w 
and the length of is at most For each Ui of ip, let Ui & V he such 

that Ui E Ui. Notice that it is possible that Ui = Ui+i for some i. Let \& be 
the subsequence of U = Ui, . . . ,Uk = W that consists in those U such that 
i < k and Ui ^ Ui^\. By definition, the sequence \i/ is a path in SjV and it 
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consists of at most k' < k nodes. The length of \l/ is at most 

fc'-l k-l 

UiUi+i\ < t\uw\ = t\UW\ 

i=l 1=1 

which means that \1/ is a t-spanning path for {U, W) in S/V. □ 



6 The Additively Weighted Delaunay Graph 

Lee and Drysdale [23j studied a variant of the Voronoi diagram called the 
Additively Weighted Voronoi diagram, which is defined as follows: Let P be 
a weighted point set. The Additively Weighted Voronoi diagram of P is a 
partition of the plane into \P\ regions such that each region contains exactly 
the points in the plane having the same closest neighbor in P according to the 
additive distance. In other words, the Voronoi cell of a pair {pi,ri) contains 
the points p such that d{p,pi) is minimum over all other pairs in P. The 
Additively Weighted Delaunay graph (AW-Delaunay graph) is defined as the 
face-dual of the Additively Weighted Voronoi diagram. 

Alternatively, if all are positive and for all i,j, we have \piPj\ > Tj + rj, 
then the pairs {pi,ri) can be seen as disks Di of radius centered at pi and 
d{p, Di) is the minimum \pq\ over all q G Di. For a set V of disks in the plane, 
we denote the AW-Delaunay graph computed from V as Del(D). When no 
two disks intersect, the AW-Delaunay graph is a natural generalization of the 
Delaunay graph of a set of points. We say that two disks A and B properly 
intersect if |y4 fl i?| > L 

Proposition 6.1 Let V he a set of disjoint disks in the plane, and A,B^V. 
The edge {A, B) is in Del(P) if and only if there is a disk C that is tangent 
to both A and B and does not properly intersect any other disk in T>. 

Proof: Suppose {A,B) is in Del(T'), and let c be a point on the boundary 
of the Voronoi cells of A and B and r be the distance from c to A. Since 
c is equidistant from A and B, it is also at distance r from B. This means 
that the disk C centered at c is tangent to both A and B. This disk cannot 
properly intersect any other disk of V, since this would contradict the fact 
that c is in the Voronoi cells of A and B. Similarly, if there is a disk that is 
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Figure 9: The Additively Weighted Delaunay graph compared with the De- 
launay graph of the disks centers. 



tangent to both A and B but does not properly intersect any other disk of 
V, then A and B are Voronoi neighbors. □ 

Note that the Additively Weighted Delaunay graph is not necessarily 
isomorphic to the Delaunay graph of the centers of the disks (see Figure |9|. 
When all radii are equal, however, the two graphs coincide. We now show 
that if P is a set of disks in the plane, then Del(T') is a spanner of "D. The 
intuition behind the proof is the following: we show the existence of a finite 
set of points P such that K{P)/V (where K{P) is the complete graph with 
vertex set P) is isomorphic to the complete graph on T> and Del(P)/P is a 



subgraph of Del(P). Then, we use Lemma 5.3 to prove that Y)e\{P)/T> is a 



spanner of which implies that Del(X') is a spanner of T>. 

Definition 6.2 Let A,B be disjoint disks and S a set of points such that 
A n S = ^ and B n S = 0. A set of points R represents S with respect 
to A and B if for every disk F that is tangent to both A and B, we have 
FnS^^^FnR^^. If V is a set of disjoint disks, then a set of points 
TZ represents V if for all A,B,C G V, there is a subset ofTZ that represents 
C with respect to A and B. 



From here to the end of the proof of Lemma 6.6 , unless stated otherwise, let 



1. y4, B be two disjoint disks in the plane having their center on the x-axis; 

2. D{y) be the disk that is tangent to both A and B and whose center 
has y-coordinate equal to y; 
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Figure 10: Illustration of the proof of Lemma 



6.3 



3. y{D) be the y-coordinate of the center of a disk D; 

4. £i,£2 be the two hues that are outer-tangent to both A and B (respec- 
tively, from below and above); 

5. yi,y2 be such that yi < y2 and D{yi) fl D{y2) ^ ^] 

6. £ be the line through the intersection points of the boundaries of D{yi) 
and D{y2) (if D{yi) and D{y2) are tangent, then i is the unique line 
that is tangent to both D{yi) and D{y2))] 

7. T(y4, B) denote the region below £2, above ii and between A and B] 
and 

8. (/~) be the closed half-plane above (below) a non- vertical line /. 

Throughout this section, it is implicitly assumed that D{oo) = ^2 ^^"^ 
D(-oo) =£^. 

Lemma 6.3 Given yi < y2 and D{yi) fl D{y2) 7^ 0, we have D{yi) fl £^ C 
D{y2) n £+ and D{y2) n r C D{yi) n r ^see Figure\l^. 

Proof: Notice that either D (?/i ) n C D (7/2 ) H or D (?/2 ) H C D (?/i ) n . 
Therefore, all we need to show is that (-D(?/2) H £^) \ (-D(yi) H £+) is not 
empty. Let ci, C2 be the respective centers of D{yi) and D{y2), and p be the 
intersection point of the infinite ray from ci through C2 with the boundary 
of D(yi)UZ}(y2). 

We show by contradiction that p is not in D{yi). If that was the case, 
then D{y2) would be completely contained in D{yi). The reason for this is 
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Figure 11: Illustration of the proof of Lemma 



6.4 



that there is no point of D{y2) that is farther from Ci than p. Let g be a 
point of D{y2). Then |gci| < \qc2\ + |c2Ci| < \pc2\ + \c2Ci\ = \pci\. But the 
fact that D{y2) is completely contained in D{yi) contradicts the fact that 
they are both tangent to A and B. 

Therefore, since p G we have p G (-0(2/2) H \ (-D(?/i) H which 
imply that D{yi) n C D{y2) n £+. Similarly, D{y2) ni^ C D{yi) n £~. □ 



Lemma 6.4 Letpi,p2 be the intersection points of the boundaries of D[yi) 
and D{y2) (if D{yi) and D {112) are tangent, then pi =P2)- Then pi and p2 
are in £2 and in if (see Figure 11). 



Proof: Let gi, q2 be the tangency points of D{yi) with A and B and Si, S2 
be the tangency points of D{y2) with A and B. By Lemma 6.3, qi,q2 are 
below i and Si, S2 are above £. Since i is above gi and q2, which are in turn 
above ii, it follows that pi and p2 are above ii. By a symmetric argument, 
Pi and are below £2- CH 



Lemma 6.5 T/ie following are true: 

1. For all p E £2 ! there exists a line y = yo such that for all disk E that 
is tangent to both A and B, if the center of E is above yo then p E E. 

2. For all p E ii , there exists a line y = yi such that for all disk E that 
is tangent to both A and B, if the center of E is below yi then p E E. 
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3. For all p in T{A, B), there exists two lines y = yo and y = yi such that 
for all disk E that is tangent to both A and B, p E E if and only if the 
center of E is between yo and yi . 

Proof: For (1), the existence of yo is guaranteed by the fact that hnij^^oo D{y) 
£2- Now, let yo be such that p G D{yo) and y' > yo- Let L{yo) and L{y') 
be the lunes respectively defined by the intersection of D{yo) and -D(y') with 



the half-plane above £2- By Lemma 6.4, the two points where the bound 



aries of D(yo) and D{y') intersect are below £2- Therefore, we have either 



L{yo) C L{y') or L{y') C L{yo). But since y' > yo, by Lemma 6.3 we have 
L{yo) C L{y') and therefore p G L{y'). The proof of (2) is symmetric. 

For (3), the existence is easy to show. Without loss of generality, assume 
d{p, A) < d{p, B). Let D be the disk centered at p that is tangent to A and 



let q be the tangency point of A and D see Figure 12 Since q G T[A, B 
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Figure 14: The five regions for Lemma 



6.6 



tliere exists y sucli tliat D{y) r\ A = q. Since D C D{y), tliere exists a disk 
tliat is tangent to both A and B and contains p. 

We now show that yi < y2 < ys implies D{yi) fl D{y^) C D{y2) (see Fig- 
ure [13]). Let £3 be the line through the intersection points of the boundaries 
of D{yi) and D{y2) and let £4 be the line through the intersection points of 
the boundaries of -D(?/2) and D^y^). Let p G D{yi) nD^y^). Since £4 is above 
£3 in D{yi) fl D^y^), p is either above £3, below £4 or both. If p G i^, then 
since ?/i < y2, by Lemma 6.3 we have that D{yi) n C D{y2) H £3 a nd 



p G D{yi) r\D{y2). Similarly, if p G ^4 , then since 2/2 < Z/s, by Lemma 6.3 



we 



have that D{y^) fl £4 C D{y2) fl £4 and p G D{y3) fl D{y2). In either case, 
p G D{y2), which completes the proof. □ 



Lemma 6.6 Let C be a disk that is disjoint of both A and B. There exists 
a set of at most six points that represents C with respect to A and B. 



Proof: Let 

C2 
C3 

C5 



{c n it) \ 

CnT{A,B) 

{Cn£tni^)\T{A,B) 



These five regions partition the disk C (see Figure 14). We show that for 



each region, there is a finite set of points that represents it. The cardinality 
of the union of the sets is no more than six. 
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Figure 15: Case C4 of the proof of Lemma 



6.6 



If Ci 7^ 0, then let yo be the minimum y such that D{y) intersects Ci. 
Let pi & Ci n D{yo). By definition of yo, for any disk E that is tangent to 
both A and B and intersects Ci, we have y{E) > yo, and by Lemma 6.5, we 
have pi E E. 

Similarly, if C2 7^ 0, then let yi be the maximum y such that D{y) 
intersects C2. Let p2 G C2 H D{yi). By definition of yi, for any disk E that 
is tangent to both A and B and intersects C2, we have y{E) < yi, and by 



Lemma [6. 5| we have p2 & E. 

If C3 7^ 0, then let yo be the minimum y > such that -D(i/) intersects C3 
and j/i as the maximum y < such that -D(?/) intersects C3. Let G C3 fl 
D{yo) and p4 G CsflZ^d/i). By definition of yo, for any disk E with > 
that is tangent to both A and B and intersects C3, we have > yo, and 
by Lemma 6^ we have ^3 G E. The same reasoning applies to p4 when 
?/(^) < 0. 

If C4 7^ 0, then let yo be the minimum y such that -D(?/) intersects C4 
and ?/i as the maximum y such that D{y) intersects C4. Let G C4 fl D{yo) 
and p6 G C4 n D{yi). Let ?/* be such that C C D{y*) (see Figure [Isj). Let 
be a disk that is tangent to both A and B and intersects C4. We show 
that y{E) < y* =^ Pb & E (and similarly, y{E) > y* =^ Pq G E). It 
is sufficient to show that y" < y' < y* ^ CD D{y") C C n D{y'). Let 
p G D{y") n C. By Lemma 6.5, 3yo{p),yi{p) such that V disk E tangent 
to both A and B, we have yo{p) < !/(-£') < l/i(p) P G -E. Therefore, the 
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Figure 16: Proof of Lemma 



6.8 



following hold: 

yo{p)<y* < yi{p) 
yo{p) < y" < yi{p) 

But since y" < y' < y*,we have y" < y' < y*, which imply that p G CnD{y'). 

Finally, since C5 fl = for any disk E that is tangent to both A and 
B, there is no need to select representative points for C5. □ 



Careful analysis of the proof of Lemma |6.6| allows us to observe that 
in fact, only two points are necessary to represent a disk C with respect 
to two other disks A and B. First, note that C4 7^ =^ C3 = and 
C3 7^ =^ C4 = 0. This reduces to four the number of points that are 
necessary. Also, if Ci 7^ and C4 7^ 0, then is on £2 and is not required 
since any disk that contains it also intersects Ci and therefore contains pi. 
Similarly, if C2 7^ and C4 7^ 0, then p^ is not required since any disk 
that contains it also intersects C2 and therefore contains p2- Therefore, if 
C4 7^ 0, then the number of points that are necessary is at most two. A 
similar argument applies to the case where C3 7^ 0. Finally, if both C3 and 
C4 are empty, then only pi and p2 may be required. Therefore, we have the 
following corollary: 

Corollary 6.7 Let V be a set of n disjoint disks. There exists a set of at 
most 2(g) points that represents V. 



Lemma 6.8 Let A and B be two disjoint disks and C be a disk intersecting 
both of them. Then there exists a disk G inside C that is tangent to both A 
and B. 
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Proof: We show how to construct G. Let a,b,c and rA,rB,rc respectively 
be the centers and radii of A, B and C. Without loss of generality, assume 
\ac\ — rc < \bc\ — tb- Let F be the disk centered at c and having radius 
rp = \bc\ — Tb (see Figure 16). The disk F is tangent to B. If F is also 
tangent to A, then let G = F and we are done. Otherwise, F is properly 
intersecting A. In that case, let p be the tangency point of F and B, I be 
the line through b and c, and G be the disk through p having its center on 
/ and tangent to A. The result follows from the fact that G is tangent to B 
and inside G. □ 



Definition 6.9 Let A and B be two disks in the plane. The distance points 
of A and B are the two ends of the shortest line segment between A and B 
(see Figure 11). If V is a set of disjoint disks, then the set of distance points 
ofV is the set containing the distance points of every pair of disks in T>. 



Theorem 6.10 Let V be a set of n disjoint disks. Then Del(I') is a t- 
spanner ofT>, where t is the spanning ratio of the Delaunay triangulation of 
a set of points. 
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Figure 19: Even if the embedding of the AW-Delaunay graph that consists 
of straight hne segments between the centers of the disks is not necessarily a 
plane graph, it is planar. 



Proof: By Corollary 6.7, let P be a set of size at most 2(g) that repre- 



sents V, let S be the set of distance points of V, and let P = RU S. Since 
Del(P) is a t-spanner of P, by Lemma 5.3, we have Del(P)/P is a t-spanner 
of K{P)/V, where K{P) is the complete graph with vertex set P. Since 
P contains the distance points of V, K{P)/T> is isomorphic to the com- 
plete graph defined on V. We show that each edge (A, B) of Del(P)/r' is in 
Y)e\{T>). Let {A,B) be an edge of Del(P)/P. This means that in P, there 
are two points a and h with a G A, 6 G P such that there is an empty circle 
C through a and h. By Lemma 6^, C contains a disk G that is tangent to 
both A and B. The disk G is a witness of the presence of the edge (A, B) 
in Del(P). If that was not the case, this would mean that there exists a disk 
F such that G n P 7^ 0. By definition of P, this implies that G n P 7^ 
and thus C (1 P ^ ^, which contradicts the fact that C is an empty circle. 
Therefore, the edge {A,B) is in Del(P). Since Del(P)/X' is a t-spanner of 
V and a subgraph of Del(X'), we conclude that Del(X') is a t-spanner of V. □ 



Note that the embedding of the AW-Delaunay graph that consists of 
straight line segments between the centers of the disks is not necessarily a 
plane graph (see Figure 19). However, the Voronoi diagram of a set of disks 
V, denoted Vor(X>), is planar [26] . Since Del(X>) is the face-dual of Vor(r'), it 
is also planar. An important characteristic of the Delaunay graph of a set of 
points regarded as a spanner is that it is a plane graph. Therefore, a natural 
question is whether Del(T') has a plane embedding that is also a spanner. 

The proof of Theorem 6.10 suggests the existence of an algorithm allowing 
to compute such an embedding: compute the Delaunay triangulation of the 
set P that contains the distance points and the representative of T>. The 
graph Del(P) can be regarded as a multigraph whose vertex set is V. Then, 
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for each pair of disks that share one or more edges, just keep the shortest of 
those edges. This simple algorithm allows to compute a plane embedding of 
Del(P) that is also a spanner of V. However, its running time is 0{n^ logn). 
Whether or not it is possible to compute a plane embedding of Del(P) that 
is also a spanner of P in a better running time remains a open question. 

7 Conclusion 

In this paper, we showed how, given a weighted point set where weights are 
positive and \piPj\ > Tj + rj for all i ^ j, it is possible to compute a (1 + e)- 
spanner of that point set that has a linear number of edges. We also showed 
that the Additively Weighted Delaunay graph is a t-spanner of an additively 
weighted point set in the same case. The constant t is the same as for the 
Delaunay triangulation of a point set (the best current value is 2.42 [T9]). 
We could not see how the Well-Separated Pair Decomposition (WSPD) can 
be adapted to solve that problem. The first difficulty resides in the fact that 
it is not even clear that, given a weighted point set, a WSPD of that point 
set always exists. Other obvious open questions are whether our results still 
hold when some weights are negative or \piPj\ < ri + rj for some i j- Also, 
we did not verify whether our variant of the Yao graph can be computed in 
time 0(?T,log?T,). Finally, another problem that could be explored is whether 
it is possible to compute t-spanners for multiplicatively weighted point sets. 
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